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Abstract. We consider the prospects for the indirect detection of low mass dark matter which couples dominantly to quarks. If the 
center of mass energy is below about 280 MeV, the kinematically allowed final states will be dominated by photons and neutral 
pions, producing striking signatures at gamma ray telescopes. In fact, an array of new instruments have been proposed, which would 
greatly improve sensitivity to photons in this energy range. We find that planned instruments can improve on current sensitivity to 
dark matter models of this type by up to a few orders of magnitude. 


INTRODUCTION 

There are three reasons to focus attention on indirect detection signatures in the photon spectrum arising from the 
annihilation or decay of very low-mass dark matter (mx ~ (9(100 MeV)) [1]: 

• If dark matter couples dominantly to quarks and is sufficiently light, then the annihilation/decay products are 
dominated by ys and 7r°s, producing striking signatures in the photon spectrum; 

• A variety of new astronomical instruments which can vastly improve sensitivity to (9(1 - 100) MeV photons are 
under consideration; 

• Direct detection experiments tend to lose sensitivity at such low masses, while indirect detection strategies may 
well be the most effective. 

For a dark matter initial state with center-of-mass energy sfs < 2m„±, the only kinematically accessible non- 
leptonic two-body final states are yy, yjfi and Since the tt® will decay almost exclusively to yy, low-mass dark 
matter is essentially “guaranteed” to produce interesting indirect detection signatures in the photon spectrum, provided 
it couples to quarks. Moreover, since sensitivity to photons in this energy range is relatively poor, there are a variety 
of proposals for new instruments which would dramatically improve sensitivity in this “MeV gap,” [2] including 
the Advanced Compton Telescope (ACT) [3], the Advanced Pair Telescope (APT) [4], the Gamma-Ray Imaging, 
Polarimetry and Spectroscopy (GRIPS) detector [2], the Advanced Energetic Pair Telescope (AdEPT) [5], the Pair- 
Production Gamma-Ray Unit (PANGU) [6], the Compton Spectrometer and Imager (COSI) [7], and ASTROGAM [8]. 
In this proceedings contribution, we describe the bounds on dark matter annihilation or decay to these final states 
arising from current data (from COMPTEL, EGRET and Eermi), and the future sensitivity of other instruments under 
consideration. 

We begin by describing a class of low-mass dark matter models which produce distinctive photon spectra. We 
then describe the general features of indirect detection strategies for probing these models, utilizing either the diffuse 
gamma ray flux or a search of dwarf galaxies. We describe the current bounds arising from existing gamma ray data 
sets in this energy range (as well as from the Cosmic Microwave Background (CMB)), and the sensitivity of future 
contemplated instruments. We end with a discussion of our results. 

MODELS AND SPECTRA 

We consider the case of low-mass dark matter whose dominant Standard Model coupling is to quarks. Indirect detec¬ 
tion is possible through either dark matter annihilation or decay; in either case, we refer to the center-of-mass energy 



of the relevant process as ^fs, and assume ^fs < 2m„±. The only non-leptonic two-body final states that are kinemat¬ 
ically allowed are yy, yjf^, and Additional leptonic or three-body final states are also available, but the relevant 
annihilation/decay rates are suppressed by factors of a or sGf, and we assume they are negligible. 

For each two-body final state, the energies of the prompt photons and neutral pions are determined by ^fs. 
The neutral pion decays through the process jf’ —> yy with ~ 99% efficiency, resulting in a “box-like” spectrum, 
characteristic of monoenergetic photons from an isotropically boosted source. As a result, all of the final state products 
produce distinctive photon signatures. 

If we assume that weak interactions are negligible and that the fundamental dark matter-quark interaction is C- 
invariant, then we may classify final states by their C eigenvalues. The photon injection spectrum associated with each 
final state is determined by kinematics, and are as follows: 

• 7r®7r° (C-even): 


• yTT® (C-odd): 


• yy (C-even): 
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Depending on the details of the model, as well as the C eigenvalue and energy of the initial state, it is possible for 
any final state to dominate. We will therefore analyze the detection prospects for each final state channel assuming the 
branching fraction to that channel is 1. 

The above photon injection spectra must be convolved with an energy resolution function to determine the spec¬ 
trum observed at any particular detector. For simplicity, we will assume in our analysis that the energy resolution 
function for each detector can be well-approximated by a Gaussian with an energy-independent fractional width. 

For channels in which a is produced just above threshold, the tt® will decay nearly at rest, producing a very 
narrow box-like feature (AE small) which may mimic a line. This feature will provide an excellent signature for a 
detector whose energy range is large enough to contain the peak. In the yV’ channel, this box-like feature will dominate 
the sensitivity of an experiment only if it is sufficiently narrow; for large enough the monoenergetic photon will 
be a larger driver of sensitivity. 


SEARCHES OF THE DIFFUSE FLUX AND OF DWARF GALAXIES 

We will focus on signals of dark matter annihilation/decay in the diffuse gamma ray flux, and in gamma ray emission 
from dwarf spheroidal galaxies. This analysis is described in detail in [1]. The diffuse flux of gamma rays has been 
measured in the energy range of interest by COMPTEL and EGRET [9]. The data are relatively smooth and can be fit 
to a power law [1]: 


d^^ldEdO. = (2.74 X 10^^ cm^2g-igj-i MeV^')(E/MeV)^^'’ 


(4) 



A indirect detection search in the diffuse gamma ray flux is essentially a search for a sharp feature in the observed 
photon spectrum. On the other hand, the observed diffuse photon spectrum constitutes the foreground to a search for 
dark matter in dwarf galaxies. 

To determine the sensitivity of current and future instruments to searches for diffuse gamma ray emission, we 
consider both a “conservative” and “optimistic” analysis. In the conservative analysis, we exclude models for which 
the expected number of photons from dark matter annihilation/decay exceeds the observed number of counts by 2cr 
in any one energy bin.' This analysis would be appropriate if there is little confidence in one’s understanding of the 
underlying astrophysical backgrounds. But in a more optimistic case, one might be confident that the underlying as- 
trophysical background exhibited no sharp features. In that case, one could instead estimate the systematic uncertainty 
in a smooth fit to the diffuse gamma ray flux (following Reference [9], we estimate the systematic uncertainty to be 
15%), and exclude models for which the expected number of photons from the dark matter signal exceeded twice this 
uncertainty in any one energy bin. 

It is important to note that, for either a conservative or optimistic analysis of the diffuse gamma ray flux, statistical 
uncertainties are subleading. As a result, sensitivity is not significantly affected by the exposure or angular resolution 
of the instrument. Instead, the sensitivity to a line signal is proportional to e *, which determines the instrument’s 
ability to distinguish a sharp feature in the data. 

For an analysis of dwarf galaxies, however, we assume that it is possible to estimate the foreground photon 
spectrum arising from astrophysical sources (and from diffuse gamma ray emission due dark matter annihilation/decay 
along the line of sight) by looking slightly off-axis. One then excludes models for which the number of photons 
expected from dark matter annihilation/decay within the dwarf galaxy could not be accommodated by a 2cr downward 
statistical fluctuation in the number of foreground photons within the set of energy bins which encompass the entire 
dwarf galaxy dark matter signal. Since the sensitivity is controlled by statistical uncertainties, it is proportional to 
•^AeffTobs/f for a line search, where is the effective area and is the run time. 

We thus see that a larger Ag^ and larger would improve the sensitivity of a dwarf galaxy search, while an 
improved energy resolution would improve the sensitivity of both dwarf galaxy and diffuse emission searches. The 
energy range, energy and angular resolutions (Icr), and effective area of the relevant instruments are listed in Table 1. 
Note that the listed values for the resolutions and effective areas are given at a particular energy, and our analysis 
is performed with these benchmark numbers. If possible, we use values corresponding to photon energies at or near 
lOOMeV with photon detection at normal incidence. For the experiments that cover a lower energy range, we use the 
performance goal values. In the case that only an estimated range is provided, we use the more pessimistic (worse 
resolution, smaller effective area) end of that range. We list two benchmark points for ASTROGAM: one for photon 
detection via Compton scattering (below lOMeV) and one for detection via pair production (above lOMeV). 


TABLE 1. Experimental parameters used to determine indirect detection bounds for dark matter annihilation 
and decay. The values between experiments are not directly comparable, since they may have been taken at 
different photon energies. The width of the point spread function (PSF) and Agjp values are not needed in the 
analysis for COMPTEL, EGRET, and Eermi, but they are included for completeness. 
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GRIPS 
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Eermi-LAT 

[13] 
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2° 
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GAMMA-400 

[14] 

100-3 ■ 10® 
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2° 

3000 

ASTROGAM (below lOMeV) 

[8] 

0.3-10 

1% 

0.5° 

119 

ASTROGAM (above lOMeV) 

[8] 

10-3000 

30% 

0.5° 

514 


’ In determining bounds from current data, we utilize the actual flux measurements and energy binning of the relevant instrument. In estimating 
the sensitivity of future experiments, we assume that the observed flux is consistent with Equation 4, and adopt an optimal energy binning for the 
relevant experiment, given an estimate for its energy resolution. 



CURRENT BOUNDS AND FUTURE SENSITIVITY 


In Figure 1, we plot current constraints on diffuse emission of photons from dark matter annihilation (left panel) and 
decay (right panel) using data from COMPTEL, EGRET and Eermi-LAT, applying a conservative analysis (the details 
are described in [1]). The yy, yn^ and channels are depicted with solid, dashed and dotted contours, respectively. 
Note that, although Eermi’s energy range goes down as far as 20MeV, their diffuse flux analysis only extends down 
to lOOMeV, in order to reduce backgrounds from cosmic ray interactions with the atmosphere. The /-factors used in 
estimating the dark matter signal flux are obtained from [15]. We also plot bounds on the annihilation cross section 
(assuming i-wave annihilation) and decay rate [16] which arise from measurements by Planck [17] of the Cosmic 
Microwave Background (CMB) spectrum. Note that, although CMB constraints exceed those of indirect detection for 
the case of dark matter annihilation, the CMB constraints are subleading for the case of dark matter decay. This is not 
surprising, given that the dark matter annihilation scales quadratically with dark matter density; the rate of dark matter 
annihilation in the early Universe is thus greatly enhanced relative to the current epoch. 


Annihilation Constraints 




FIGURE 1. Diffuse gamma ray flux constraints on dark matter annihilation (left panel) and decay (right panel) to the yy (solid), 
yTfi (dashed), and (dotted) channels obtained from data from COMPTEL (blue), EGRET (red) and Eermi-LAT (green). Also 
shown are constraints which can be obtained from Planck data (black). The vertical black dashed lines are the kinematic thresholds 
for the yrr** and channels. 


In Eigure 2, we plot the sensitivity of a variety of future instruments, in the case of either a diffuse emission search, 
or a 5-year search of photons from Draco (the /-factor from Draco is obtained from [18]). A detailed description of 
the analysis is found in [1]. The results of conservative and optimistic analyses of diffuse emission are depicted as the 
boundaries of the solid bands, while the sensitivity of the analysis of Draco is depicted by the hatched regions whose 
thicknesses represent the Icr systematic uncertainties in /. 

One can see several features in the sensitivities shown in Eigure 2. In the yy and yuP channels, ASTROGAM’s 
Compton detector can detect the lower energy photons with a high energy resolution; consequently, it has the greatest 
sensitivity in the yy channel in a diffuse gamma ray search. Similarly, the enhanced sensitivity of GRIPS over AdEPT 
(and ASTROGAM in the pair-production regime) to the yy channel in a diffuse gamma ray search is due to its better 
energy resolution. ASTROGAM is the least sensitive for a diffuse search in its pair production regime, given the 
pessimistic projected 30% energy resolution; using a more optimistic 20% would bring its sensitivities more in line 
with those of AdEPT. However, AdEPT, GRIPS, and ASTROGAM all have similar sensitivity to the same channel 
in a search of dwarf galaxies because the greater effective area of AdEPT and ASTROGAM compensates for poorer 
energy resolution. Eor the yy and yjfi channels, the sensitivity of GRIPS is suppressed at higher energies due to the 
the relatively low upper limit of its energy range (80MeV), compared to AdEPT and ASTROGAM. Additionally, 
GRIPS will miss a fraction of the box-like feature arising from decay, unless the pion is nearly at rest. 

Comparing the cases of dark matter annihilation and dark matter decay, one sees that a search of dwarf spheroidal 
galaxies shows a greater improvement in sensitivity, relative to diffuse flux searches, for the case of dark matter 
annihilation. This is not surprising, since the enhanced dark matter density within a dwarf galaxy will lead to a 
quadratic enhancement in the annihilation rate, but only a linear enhancement in the decay rate. 






















CONCLUSIONS 


We have considered a class of models in which low-mass ( < 2m„±) dark matter couples dominantly to quarks. For 

these models, kinematic considerations favor the 77, 77r® and tt'V' final states for dark matter annihilation or decay, 
regardless of the details of the interaction of dark matter with quarks. These models thus generically predict sharp and 
striking signals in the photons spectrum in the energy range (9(1 - 100) MeV. 

This is especially interesting because, although there is currently a dearth of instruments focused on this en¬ 
ergy range, the astrophysics community is considering several plans to fill this “MeV gap.” These new instruments 
could improve sensitivity to these indirect detection channels by up to a few orders of magnitude. Indeed, although 
constraints from CMB observations are by far the most stringent current limits on dark matter annihilation, searches 
of dwarf galaxies by future instruments could exceed those limits by more than an order of magnitude. New instru¬ 
ments searching for photons in the (9(1 - 100) MeV range may exceed previous instruments in effective area, energy 
resolution and angular resolution, providing a variety of paths for improving sensitivity to low-mass dark matter. 
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FIGURE 2. The projected sensitivity of ASTROGAM, AdEPT, and GRIPS to dark matter annihilation (left panels) or decay (right 
panels) in the 77 (top row), 721 ® (middle row) and (bottom row) channels. The shaded regions are bounded by contours which 
describe sensitivity of a conservative or optimistic analysis of diffuse gamma ray emission, as described in the text. The hatched 
region describes the sensitivity contour for a search of Draco, when the /-factor is varied through Icr systematic uncertainties. The 
black line describes the exclusion contour from Planck data. The dashed vertical lines are the kinematic thresholds for the y/r® and 
;r®;r® channels. 






























































